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Abstract 



We present data on multiple production of jets with transverse energies near 
20 GeV in pp collisions at -s/s = 1.8 TeV. QCD calculations in the parton- 
shower approximation of PYTHIA and HERWIG and the next-to-leading order 
approximation of .TETRAD are compared to the data for one, two, three, and 
four jet inclusive production. Transverse energy spectra and multiple jet angu- 
lar and summed transverse-energy distributions are adequately described by 
the shower approximation while next-to-leading order calculations describe 
the data poorly. 
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I. INTRODUCTION 



The study of multiple jet production at high transverse energy was a goal of the 1993- 
1995 run of the Fermilab Tevatron collider, and the results have been compared with leading- 
order QCD predictions by both the CDF [1] and D0 [2] collaborations. These high-E'^ data, 
where Et is the transverse energy of the jet, are described satisfactorily by complete tree- 
level leading order 2 ^ N QCD calculations [3] and by the HERWIG parton-shower Monte 
Carlo [4] program. This kinematic region is described by Q^/s ~ 1, where is the square of 
the momentum transfer between partons (which we set equal to £"1^), and s is the square of 
the partonic center of mass energy. In this paper, we describe jet production measurements 
at significantly lower values of Et where detailed measurement of jet production in this 
kinematic region can provide information on the evolution of higher-order jet processes. 
In the same low Et region the D0 collaboration has previously reported the ratio of the 
inclusive three-jet to the inclusive two-jet cross section as a function of the scalar sum 
of jet transverse energies {Ht — Y.Et) with Et > 20 GeV [5]. The ratio data can be 
described by the jetrad next-to- leading order Monte Carlo [6] program. In this paper we 
make comparisons between Monte Carlo and several characteristics of multiple jet events 
including the leading jet transverse energy, the relative azimuthal angle between jets, and 
the summed vector transverse momenta of jets. 

II. DATA SAMPLE AND CORRECTIONS 

The data were collected with the D0 detector at a proton-antiproton center-of-mass 
energy ^/s = 1.8 TeV. Jets were identified using the liquid- argon uranium calorimeters, 
which have segmentation of Ai] x A</) = 0.1 x 0.1, where pseudorapidity r] = — lntan|, 
6 is the polar angle, and (p is azimuthal angle [7]. At least one calorimeter trigger tower 
(A77 X A0 = 0.2 X 0.2) with Et > 2 GeV was required by a hardware trigger, and at least 
one jet with Et > 12 GeV was required by a subsequent software trigger [8]. Jets were 
reconstructed using a fixed cone algorithm with radius ATZ = \/ Ar]'^ + = 0.7 in 77 — 
space [8]. The jet reconstruction threshold was Et = 8 GeV. If two jets overlapped and the 
shared transverse energy was more than 50% of the transverse energy of the lower- i^^ jet, 
the jets were merged; otherwise they were split into two jets. The integrated luminosity of 
this data sample is 2.0±0.3 nb~^. Instantaneous luminosity was restricted to be below 3 
X 10'^° cm^^s^^ to minimize the number of multiple pp interactions in a single beam crossing. 

To provide events of high quality, online and offline selection criteria suppressed multiple 
interactions, the cosmic ray background, and spurious jets [8]. Jets were restricted to the 
pseudorapidity interval |?7| < 3. The primary vertex of each event (reconstructed from time- 
of-flight as measured by scintillation counters [7]) was required to be within 50 cm of the 
detector center. 

Jet energies have been corrected for calorimeter response, shower development, and 
various sources of noise [9]. These corrections constitute the largest source of systematic 
uncertainty on the jet cross section. Typical values of the jet energy correction are (15- 
30)%, with an uncertainty of (2-4)%. In our study, we consider jets with Et > 20 GeV; 
for an inclusive n-jet event, the n jets with the maximum Et (the leading jets) must have 



5 



transverse energy above the threshold value. For example, a 3-jet event must have at least 
3 jets above 20 GeV. The trigger efficiency is 0.85 for the inclusive (n = 1) jet sample for 
energies near threshold, rising rapidly to unity at larger Et- The efficiency is essentially 
unity for n > 1. 

To compare with data, Monte Carlo (MC) events were generated using the pythia 6.127 
[10], HERWIG 5.9 [4], and jetrad [6] programs, pythia and herwig simulate particle-level 
jets in the parton-shower approximation. JETRAD simulates jets in the next-to-leading order 
approximation. To simulate detector resolution effects, the MC jet transverse energies were 
smeared with the experimentally determined jet energy resolution [9], which is ~ 20% at 
Et — 20 GeV. Jet angular smearing used rj and resolutions obtained by a MC simulation 
of the calorimeter response using HERWiG 5.9 and GEANT [11]. These resolutions are ~ 0.08 
at Et = 20 GeV. In pythia and herwig, jets were reconstructed at the particle level using 
the D0 algorithm, and in jetrad, at the parton level, using the Snowmass algorithm [12]. 

III. LEADING JET Et DISTRIBUTIONS AND SYSTEMATIC UNCERTAINTIES 

Distributions in transverse energy for the leading jet for inclusive n—1 to n—A jet events 
are shown in Fig. 1, along with the results from pythia and HERWiG simulations. In these 
and all other plots, the data have been corrected for inefficiencies and energy calibration, 
but not for contributions from the underlying event. All simulated distributions have been 
smeared with energy and angular resolutions. Also to describe the data quantitatively, we 
normalize the theory (with a factor of 0.75 for pythia and 1.6 for herwig) to the observed 
two-jet inclusive cross section in Fig. 1(b) for Et > 40 GeV. 

The normalised theory is in agreement with the data for all of the jet samples over the en- 
tire Et interval. A detailed comparison is shown in Figs. 2 and 3. Here the simulations have 
been brought into agreement with the data by selecting parameters that enhance low Et 
jet production. In the case of pythia, the core of the hadronic matter distribution [10] has 
been increased to the fraction 0.32. An increased core fraction (the parameter PARP(83)) 
leads to enhancement of the multiple interaction rate [10], which tends to produce events 
with large multiplicity because of additional radiated low energy jets and underlying event 
energy. In the case of herwig, the minimum transverse momentum for the hard subpro- 
cesses has been set to 3.7 GeV. A decreased minimum transverse momentum (the parameter 
PTMIN) leads to increased soft underlying event contributions. The default values for these 
parameters are PARP(83)=0.5 and PTMIN=10 GeV. Variation of these values by more 
than 15% leads to disagreement with the low Et data. Other parameters, when varied from 
their default values, do not change the distributions significantly. 

Figures 2 and 3 show the fractional difference (Data — MC) / MC for the Et spectra 
in Fig. 1 with the uncertainties arising from jet-energy calibration and resolutions. The 
systematic uncertainty on the cross section is due primarily to the uncertainty in the energy 
calibration. This uncertainty can be estimated by considering cross sections derived with 
±1 standard-deviation corrections to the jet energy scale. The same procedure can be used 
to derive the uncertainties due to jet energy and angular resolutions in the MC. At Et — 25 
GeV, the uncertainty in the three-jet cross section due to cahbration of the data is 39%, and 
uncertainties in the MC due to energy and angular resolutions are 19% and 7%, respectively. 
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FIG. 1. The transverse energy distributions of the leading jet for (a) single-inclusive, (b) two-jet 
inclusive, (c) three-jet inclusive, and (d) four-jet inclusive events. Solid histograms show the PYTHIA 
simulation normalized (with a factor of 0.75) to the inclusive two-jet sample for Et > 40 GeV. 
Dotted histograms are similarly normalized HEKWIG results (increased by a factor of 1.6). 
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FIG. 2. (Data — pythia)/pythia as a function of the transverse energy of the leading jet for 
(a) single-jet inclusive, (b) two-jet inclusive, (c) three-jet inclusive, and (d) four-jet inclusive event 
samples. The relative systematic uncertainties in the cross section corresponding to the energy 
calibration added in quadrature with 15% uncertainty in luminosity are shown by the solid lines. 
The uncertainty in the ratio (Data — MC) / MC from energy and angle smearing is shown by the 
dashed lines. The total uncertainty on the ratio is shown by the dotted lines. 
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FIG. 3. (Data — herwig)/herwig as a function of the transverse energy of the leading jet for 
(a) single-jet inclusive, (b) two-jet inclusive, (c) three-jet inclusive, and (d) four-jet inclusive event 
samples. The relative systematic uncertainties in the cross section corresponding to the energy 
calibration added in quadrature with 15% uncertainty in luminosity are shown by the solid lines. 
The uncertainty in the ratio (Data — MC) / MC from energy and angle smearing is shown by the 
dashed lines. The total uncertainty on the ratio is shown by the dotted lines. 
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The uncertainty from energy resohition represents the dominant uncertainty in the MC. In 
Figs. 2 and 3, the relative systematic uncertainties in the cross section corresponding to 
the energy cahbration added in quadrature with 15% uncertainty in luminosity are shown 
by the sohd hnes. The uncertainty in the ratio (Data — MC) / MC from energy and angle 
smearing is shown by the dashed lines. The total uncertainty on the ratio is shown by the 
dotted lines. As shown in Figs. 2 and 3, both pythia and HERWIG describe the data quite 
well. 

IV. TRANSVERSE ENERGY AND AZIMUTHAL DISTRIBUTIONS 

To explore features of three- and four-jet production, we turn to observations of relative 
azimuthal distributions, distributions in summed transverse momenta, and three-jet studies. 
In Fig. 4(a) we plot the azimuthal difference between the leading two jets in events with two 
or more jets. Figures 4(b-d) show the azimuthal difference between the first and second, 
first and third, and second and third highest-i?T jets in three-jet events. In Fig. 4(a) we 
see the strong anticorrelation (in the transverse plane) expected of two-jet events. The 
peak of the distribution widens substantially in the three-jet sample (Fig. 4(b-d)). The 
peaks correspond to the kinematic constraint of transverse momentum conservation for jets 
produced in hard QCD subprocesses. PYTHIA (normalized as in Fig. 1) approximates the 
observed three-jet cross section and shapes. However, small discrepancies with HERWIG (also 
normalized as in Fig. 1) are evident. 

Distributions of the square of the summed vector transverse momenta of jets Q\ = 
(Eti + H — ■ + ^TnY in Fig. 5 show significant imbalance of the transverse momenta for 
n leading jets. If events at large are removed by requiring balanced transverse energy, 
the corresponding three- and four-jet cross sections of Fig. 1 decrease at small Et- The 
shoulder at ~ 1600 GeV^ in Fig. 5(a) can be eliminated by restricting the event sample 
to just two jets with Et above 20 GeV, and no other jets between 8 and 20 GeV. This 
shoulder can consequently be associated with higher-order radiation. 

To find the pair of jets {i,j} most likely to originate from the hard interaction (rather 
than from gluon bremsshtrahlung) , we define the scaled summed dijet vector transverse 
momentum: qij = (Eti + Etj) / {ETi + Etj) ■ We choose the pair with the smallest magnitude 
of this vector and in Figs. 6(a) and 7(a) plot the distribution of the relative azimuthal angle 
$c between the jets in that pair. The data, PYTHIA, and HERWIG show a narrow maximum 
in the region where two jets from the hard scatter appear back-to-back ($c = ^r). The 
prediction from JETRAD is peaked away from <I>c ~ because only one extra jet is present. 

Figures 6(b,c) and Fig. 7(b,c) show the azimuthal separation of the third jet from each 
of the two jets that correspond to the minimum qfj. These distributions contain events only 
for TT — $c < 0.4; that is, events in which the balanced jets are essentially back-to-back. If 
the third jet were correlated with the balanced jets, it would be observed nearby or opposite 
the balanced jets. However, the data show the third jet to be weakly correlated with the 
balanced jets, and emitted at all angles. The uncertainties associated with energy calibration 
and luminosity are shown by the solid lines in Figs. 6 and 7. Uncertainties from the energy 
resolution are shown by dashed lines in Fig. 6. 
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FIG. 4. Distributions of the relative azimuthal angle between two jets in (a) two-jet inclusive 

events and in three-jet inclusive events (b-d). Jets arc ordered by their transverse energies. The 
PYTHIA predictions are indicated by the solid histograms and the herwig predictions by the dotted 
histograms. 



11 




FIG. 5. Distributions of the square of the summed vector transverse momenta Q\, for (a) 
two-jet inclusive, (b) three-jet inclusive, and (c) four-jet inclusive event samples. The PYTHIA 
predictions are indicated by the solid histograms and the HERWIG predictions by the dotted his- 
tograms. 
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FIG. 6. Azimuthal distributions between the leading jets in 3-jet events. The data is shown 
by the closed circles. Panel (a) shows the azimuthal separation between the two jets with the 
minimum summed transverse energy. Panel (b) shows the azimuthal separation between the third 
leading jet and the first jet of the minimum transverse energy pair. Panel (c) shows the azimuthal 
separation between the third leading jet and the second jet of the pair. PYTHIA is given by the solid 
histograms, JETRAD is shown by the dotted histograms. The uncertainties associated with energy 
calibration and luminosity are shown by the solid lines. Uncertainties from the energy resolution 
are shown by dashed lines. 
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FIG. 7. Azimuthal distributions between the leading jets in 3-jet events. The data is given 
by the closed circles (all jets) and by the closed triangles (the jets overlapped with more than 
one jet are excluded). Panel (a) shows the azimuthal separation between the two jets with the 
minimum summed transverse energy. Panel (b) shows the azimuthal separation between the third 
leading jet and the first jet of the minimum transverse energy pair. Panel (c) shows the azimuthal 
separation between the third leading jet and the second jet of the pair. HERWIG is given by the 
solid histograms (all jets), and the dotted histograms (the jets overlapped with more than one jet 
are excluded). The uncertainties associated with energy calibration and luminosity are shown by 
the solid lines. 
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We sec that the data, PYTHIA, and HERWIG have wider distributions than jetrad. 
PYTHIA describes the data quite well, while JETRAD fails. The agreement with PYTHIA has 
been achieved only with enhanced multiple parton interaction rates, herwig demonstrates 
small quahtative disagreement with the shape of the azimuthal plot of Fig. 7(b); the peak 
at 7r/2 is produced by jets reconstructed from the underlying-event energy [4] and grows 
quickly with small changes in PTMIN. Such jets are strongly overlapped with more that 
one jet. If jets overlapping two or more nearby jets are excluded, the HERWIG shape in 
Fig. 7(b) improves but the agreement shown in Fig. 7(a) worsens. (The cone algorithm 
reconstructs jets from seed towers and may therefore reconstruct jets sharing energy. The 
reconstruction algorithm then merges or splits the energy encompassed in these overlapping 
jets [8].) Elimination of these jets tends to suppress contributions from the soft underlying 
event. Soft interactions result in a wide distribution of particles throughout angular phase 
space. Jets reconstructed from these particles tend to be wider and of lower energy than 
more coUimated partonic jets. Such jets often share a significant fraction of energy with 
similar, neighbouring jets and are merged into a single jet. 

The shapes of the simulated distributions are sensitive to modeling of the multiple parton 
interactions. Tuning of the multiple interaction contribution in PYTHIA and the minimum 
generated transverse momentum in HERWiG are required for good agreement. In particular, 
simulations with smaller contributions from soft parton interactions show discrepancies with 
the data. 

V. CONCLUSIONS 

In this paper we showed comparisons between Monte Carlo and data for several charac- 
teristics of multiple jet events with a low jet-£^x threshold. These comparisons included the 
leading jet transverse energy, the relative azimuthal angle between jets, and the summed 
vector transverse momenta of jets. Our data on multiple jet production at low Et agree 
with PYTHIA and herwig. This is observed in the distributions of the transverse energy of 
the leading jets (Fig. 1), azimuthal distributions (Fig. 4), in the square of the summed vector 
transverse momenta Q"^ (Fig. 5), and in the three-jet angular distributions that suggest the 
presence of a weakly correlated jet (Figs. 6, 7). jetrad cannot adequately describe the 
angular distributions of the three leading jets in three jet events. 
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